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Abstract. The quantum theory of the damped harmonic oscillator has been a subject 
of continual investigation since the 1930s. The obstacle to quantization created by the 
dissipation of energy is usually dealt with by including a discrete set of additional 
harmonic oscillators as a reservoir. But a discrete reservoir cannot directly yield 
dynamics such as Ohmic damping (proportional to velocity) of the oscillator of interest. 
By using a continuum of oscillators as a reservoir, we canonically quantize the harmonic 
oscillator with Ohmic damping and also with general damping behaviour. The 
dynamics of a damped oscillator is determined by an arbitrary effective susceptibility 
that obeys Kramers-Kronig relations. This approach offers an alternative description 
of nano- mechanical oscillators and opto-mechanical systems. 



PACS numbers: 03.65.-w, 03.65.Yz, 03.70.+k 
1. Introduction 

Few classical dynamical systems are as simple or important as the one-dimensional 
damped harmonic oscillator: 

q + jq + cJ 2 q = /(*). (1) 

But the simplicity of this dynamical system does not survive the transition to quantum 
mechanics. The presence of dissipation in (or amplification with t — > —t) leads 
to severe difficulties with its quantization, a problem that has attracted repeated 
investigation from the 1930s to the present day (some historical information is given 
in [TJ [21 [3] ) . A brief account of the different approaches to quantizing Q will explain 
the rather straightforward, but crucial, respect in which the starting point of this paper 
differs from previous work. 

The central obstacle in quantizing the dynamics ([I]) is that a Hamiltonian is required 
that will generate the quantum time evolution (a Lagrangian will suffice for the path- 
integral approach). A Lagrangian and Hamiltonian exist that give the equation of 
motion (jl]), but they are time dependent because of the dissipation (or amplification) 
and this leads to difficulties in implementing the canonical commutation relation [IJ |2] . 
If the dissipated energy is included through extra dynamical degrees of freedom, then 
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the problem of a time-dependent Hamiltonian can be avoided. The simplest approach 
is to write a two-body Hamiltonian describing one damped and one amplified oscillator, 
with conserved total energy jl]. But the canonical variables of this system pljj are not the 
positions and momenta of the two oscillators, so again there is no straightforward way of 
imposing canonical commutation relations on each oscillator [H El EJ El El E] • In recent 
times the most popular approach is to treat a damped harmonic oscillator as a free 
oscillator coupled to a reservoir of oscillators of different frequencies. The apparently 
universal practice for investigations of the damped harmonic oscillator has been to use 
a discrete set of oscillators for the reservoirjf] The resulting form of the Hamiltonian 
is attributed to Magalinskii [9], and it is also the most popular starting point for 
attempts to describe quantum Brownian motion (with a free particle coupled to the 
reservoir) QUI [HI H21 EEH HH EH HS1 El HE] ■ This approach does not give the dynamics 
([TJ for the oscillator of interest, with a damping proportional to velocity; instead, the 
oscillator equation of motion can be written with a term J^ Q dsg(t — s)q(s), where g{t) 
is an integral kernel dependent on the coupling to the reservoir [3] . In order to describe 
damping proportional to velocity, often called Ohmic damping in view of the electrical 
application of Q, a limiting procedure g(t) — > 27<5(t) must be employed at some 
point [3] . The limit that produces Ohmic damping (or amplification) involves arbitrarily 
decreasing the frequency spacing between oscillators in the reservoir while maintaining 
reservoir oscillators of arbitrarily high frequencies; in other words, the reservoir is 
effectively regarded as a continuum, but only after the dynamical equations have been 
solved under the assumption that the reservoir is a discrete set. The impossibility of 
achieving Ohmic damping with a discrete reservoir and the emergence of Ohmic damping 
through a delicate continuum limit, after the dynamics has been solved, is discussed in 
great detail by Tatarskii [19]. In addition to the approaches just described, which seek to 
employ the standard quantization rules, there are phenomenological approaches to the 
damped harmonic oscillator where no rigorous quantization is attempted (reference (20] 
is just one example of such approaches). 

The starting point of the results presented here is a harmonic oscillator coupled 
to a reservoir, where the latter is a continuum of oscillators of all positive frequencies. 
Use of a continuum reservoir from the outset gives a much richer dynamical system 
compared to the use of a discrete reservoir. As the system we will be analyzing has an 
uncountable number of degrees of freedom, it acquires many of the properties of a field 
theory, and is qualitatively different from a countable set of coupled oscillators, even if 
the latter set is infinite. We will describe the general method for solving the dynamical 
system with a continuum reservoir and solve exactly the case of Ohmic damping and 
amplification, which will emerge from a particular choice of coupling to the reservoir. 
The quantization for general damping and Ohmic damping is treated in detail, including 
the diagonalization of the Hamiltonian and the case of thermal equilibrium. 

J There are of course countless other contexts in which use is made of a discrete reservoir. 
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The continuum reservoir appears to originate with Huttner and Barnett [2T], a 
paper that is very well known and yet whose technical innovation and importance 
have been under-appreciated. An example of what can be achieved with a continuum 
reservoir is the canonical quantization of the macroscopic Maxwell equations for 
arbitrary media obeying the Kramers-Kronig relations {221 123"] , including bi-anisotropic 
and moving media [2"H 125] . In the present paper the continuum reservoir will again 
show its power by allowing an exact treatment of the dynamics ([!]) and its canonical 
quantization. More importantly, the continuum reservoir will naturally lead to the 
characterization of a quantum damped harmonic oscillator by an arbitrary effective 
susceptibility obeying Kramers-Kronig relations. Because of this last property, the 
quantum damped harmonic oscillator will be found to have much in common with the 
quantum theory of light in macroscopic media. This offers a alternative framework for 
describing macroscopic quantum oscillators, which are now a subject of some remarkable 
experiments [261 123 1211 [29J, [30] . Rather than attempting to capture the immensely 
complicated microscopic physics, the results in this paper suggest that macroscopic 
quantum oscillators may be describable by effective susceptibilities that are to be 
experimentally measured, just as the electromagnetic properties of macroscopic media 
are captured by measured permittivities and permeabilities. 

The price to be paid for employing the continuum reservoir is the extra 
mathematical complexity compared to the discrete case. Many aspects of this 
mathematical apparatus, which provides an important and unusual addition to standard 
quantum field theory, have still not been fully explored. As well as treating the specific 
problem of the damped harmonic oscillator, the results presented here give further 
insight into the remarkably rich classical and quantum physics of a continuum reservoir. 

Section [2] gives the Lagrangian and equations of motion of an oscillator coupled to 
a continuum reservoir. In sections [3] and [4] the dynamics for a particular coupling that 
gives damping proportional to velocity is solved in detail. The system is quantized in 
section [5] and the diagonalization of the Hamiltonian is addressed. In sections [6] and [7] 
coherent-state solutions and thermal equilibrium are treated. Possible applications and 
extensions of the results are discussed in section [8l 

2. Lagrangian and dynamical equations 

We consider the dynamical system with Lagrangian 

1 1 1 r°° r ■ \ f°° 

L =2?- ^ + 2j dU \ X ' ~ u2X ') + J dua ^ X - ( 2 ) 

This describes a harmonic oscillator with displacement q and frequency Uq and a 
reservoir of oscillators with displacements X u and frequencies u G [0, oo), with the 
g-oscillator linearly coupled to the reservoir by an arbitrary coupling function a(u). We 
couple only displacements in (|2]), not velocities, and the reservoir oscillators are not 
directly coupled to each other. The variables q(t) and X^it) are functions only of time, 
so the entire system may be viewed as located at one point in space. The dynamical 
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variables in (|2]) are thus not fields in the conventional sense, but the dependence of 
X u (t) on the continuous quantity u will give this system many of the properties of a 
field theory, in sharp contrast to the case of a discrete reservoir. The total energy of the 
system described by ^ is 

1 1 1 f°° / • \ f°° 

E =2 < j 2 + 2 WoV + 2 J dU \ X " + U2X ") ~ J dU a ^ )qX »- (3) 

The system (|2]) is close to being a drastic simplification of the Huttner-Barnett 
model [21] of a dielectric coupled to the electromagnetic field: if the electromagnetic 
field is removed along with the spatial dependence of the medium and the reservoir, 
then the Huttner-Barnett model almost reduces to the system ([2]), the difference being 
that the g-oscillator would be coupled to X^ rather than X u . The relationship of 
our system to the Huttner-Barnett model will be commented on at several points, as 
some of the technical achievements of Huttner-Barnett will be closely related to results 
here. The main differences from Huttner-Barnett, in addition to the simplifications just 
described, are that (i) we do not rely solely on a retarded or advanced solution of the 
reservoir dynamics, (ii) much of our time will be spent in obtaining the exact solution 
for a specific coupling function, whereas Huttner and Barnett consider only a general 
coupling function, and (iii) our coupling term is of a different form. 

The Euler-Lagrange equations of Q are 



q + u^q — / da;a(u;)X w = 0, (4) 
Jo 

X u + w 2 X w - a(u)q = 0. (5) 
The general solution of the reservoir equation ^ can be written 

X u (t) = Aq(uS) cos cut + B (u) sin cut H / dt' q(t') sin [uj(t — t' 

u Jo 



e -o+|*-t'| 



A (u)=X u (0), B (u) = -X u (0), 

CO 



(6) 



where + is a positive infinitesimal quantity. By means of (J6J) we can impose arbitrary 
displacements X u (0) and velocities X u (0) on the reservoir at t — (these are not "initial" 
conditions because the solution ^ is valid for all t). The choice of t = is of course 
arbitrary, but no generality is lost by the form ^ if we wish to impose conditions on the 
reservoir at some finite time (the imposition of conditions in the infinite past or future 
is dealt with later in this section). The presence of the exponential in ^ is important 
for taking the Fourier transform of this general solution, and it can be understood as 
follows. Solutions of ^ can be constructed using a Green function G(t) defined by 

G + u 2 G = 5{t). (7) 

For example, the retarded G r (t) and advanced G a {t) Green functions are 

G r (t) = -9(t)smute-° +t , (8) 
to 

G a (t) = --9(-t) sinut e 0+t , (9) 
to 
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where 8{t) is the step function. The exponential factors are necessary in the Green 
functions (JsJ) and ^ for their Fourier transforms to exist, and the infinitesimal number 
+ gives the familiar pole prescriptions in the frequency domain, with G r (uj) analytic in 
the upper-half complex w-plane and G a (u) analytic in the lower-half plane. The general 
solution ^ is constructed with the difference 



G r {t) — G a {t) = —sinute 

CO 



-0+1*1 



(10) 



which is a solution of the homogeneous version of ^ (i.e. §7§ without the delta function). 
This gives the exponential factor in ^ that is required to define the Fourier transform 
ofX u (t). 

To complete the solution for the dynamics we must substitute ^ into Q and solve 
the resulting equation for q(t): 



q + u Q q 



poo f*t 2 / \ 

/ duj dt'qW^-smiuit-t'^e- ^-^ 
Jo Jo u 



duj a(u) < - [A (uu) + iB Q (uj)] exp(— iut) + c.c. > = 0. 



For most coupling functions a(u) this integro-differential equation is difficult to solve 
in the time domain. When written in the frequency domain however, (11) becomes 
an integral equation that can be solved by a systematic procedure once the coupling 
function a{u) is chosen. Instead of directly Fourier transforming (11), it is a little 
simpler to Fourier transform ^ and substitute the result into the Fourier transform of 
Q. Defining the relation 

1 f°° 

du f(u) exp(— iut) 



f(t) 



2tt 



(12) 



between the time and frequency domains, we obtain from (|6| 
X u (u/) = 7T [A (u) + LB (w)] S(u -u')+ir [A (u) - iB (u)) 5(u + u') 



U 



5(u + u')]P / d£ 



g(0 



(13) 



u* — u'" 2uj 

where the infinitesimal quantity + in ^ has given rise to principal-value terms denoted 
with a P (delta-function terms arising from + are found to cancel). The Fourier 
transforms of the basic equations Q and ^ are 



u' 2 + ul)q(u') 



du a(oj)X u (u') = 0, 



i-uo' 2 + u^X^u') - a(u)q(u') = 0. 



(14) 
(15) 



One verifies that (13) solves (|15j), and substitution of (13) in (14) yields an integral 
equation for q(co'): 



UJ 



Un 



uj- 1 



u' 



9M + ^p 



de 



g(0 



-7ra(\u'\)[A (\u'\)+isga(u')B (\u'\ 



(16) 
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Equation ( [16] ) is consistent with the relation q*(co') = q(-uj'), which holds because q(t) 
is real. Once q{oj') is found from (16), the reservoir solution X u (t) is most easily found 
by substituting q(t) into 

In deriving the frequency-domain result (16) we have indifferently commuted time 
and frequency integrations; if this commutation is not valid then the solution q(u') of 
(16) will not give solutions q(t) and X u (t) of the dynamical equations (HJ) and (J5|. It 
turns out that for a constant (i.e. frequency-independent) coupling a{u) = a, the time 
and frequency integrations do not commute, but this case of constant coupling is easily 
solved in the time domain. For (16) to be valid the coupling function must also be such 
that the integral containing a 2 {uj) converges. Needless to say, it must always be directly 
verified, regardless of the solution method employed, that the dynamical equations ^ 
and @ are satisfied. 

Once a coupling function a(u) and the t = state of the reservoir are chosen, 
(16) presents a principal- value integral equation for q(u'), also known as a singular 
integral equation. In the standard classification, (16) is an inhomogeneous singular 
integral equation of the third kind. There is a systematic and very elegant method of 
solving singular integral equations, which exploits some complex analysis of Riemann 
and Hilbert [3TJ, [32] . Integral equations such as (16) are solved by relating the equation 
to a boundary problem of the theory of analytic functions, usually called a Riemann 
problem but sometimes called a Hilbert problem [31, 32J. If the integral equation has a 
solution it can be found using any solution of the associated Riemann problem [3T| [32] . 
We refer the reader to Pipkin's text [31 J for the details, which are too lengthy to be 
described here (the classic text of Muskhelishvili [32J does not appear to treat directly 
equations on an infinite interval). 

The following additional remarks on the dynamical equations are worthwhile. Note 
that for any coupling function for which (16) is valid, we can insert any function 
q{uj'){= q*(—cu')) and solve for Aq{u') and B (co'), provided the integral involving q(u>') 
converges. In other words, we can choose the dynamics of the g-oscillator from an 
enormous class of functions q(t) whose Fourier transforms exist, and use (16) to find the 
values of the displacements and velocities of the reservoir at t — that will produce this 
dynamics q(t). It must be borne in mind that there may well be solutions q(t), X u {t) 
whose Fourier transforms do not exist; for example, there is no mathematical reason why 
the function q{t) should be bounded. Even if q(cu) exists for a solution q(t), it can only 
be found by solving (16) if the integral containing q(u) is well defined. Transformation 
to the frequency domain is an essential part of diagonalizing the Hamiltonian, which in 
turn gives the most convenient description of quantum states of the system. But there is 
no guarantee that the diagonalized Hamiltonian will give all the dynamical solutions of 
the original Hamiltonian, and moreover there are couplings for which the Hamiltonian 
cannot be diagonalized (see section [5]). 

It is already clear that our dynamical system is completely different from that 
in which the continuum reservoir is replaced by a discrete one. The analysis of 
dynamical systems consisting of a discrete number of coupled harmonic oscillators is 
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largely a matter of algebra [19]. In our case the mathematics is more challenging and 
this is because the dynamics is much richer. The possible couplings a(u) constitute 
the entire space of real functions, rather than a discrete set of numbers. Moreover, 
for some coupling functions a(u), the existence of a solution of (16) requires severe 
constraints on the functions A (u') and B (u'), whereas for other coupling functions 
there will be no such constraints [31] . The general mathematical theory also shows that 



there is a subspace of coupling functions for which (16) gives only the trivial solution 
q(t) = X u (t) = when X w (0) = X u (0) = [31], although it is by no means clear what 
this subspace is. It would be very interesting to know how the space of solutions behaves 
as a functional of coupling a(u), but that is far beyond the scope of this paper. We 
will however find a constraint on the coupling function that must be satisfied for the 
Hamiltonian to be diagonalizable. 

In addition to the form (|6]), we also require the general solution of the reservoir 
equation ^ in a form that allows the imposition of conditions in the infinite past or 
future. We choose the infinite past t — > — oo, and the required form of the general 
solution of ([5]) is obtained by using the retarded Green function (|8]), yielding 

a(oo) rt 



X w (t) = A R (oo) cos ut + B R (oo) smut + 



00 



dt' q(t') sin [co(t — t')} e 



-o+(t-t') 



(17) 



Here the functions A R (u) and B R (u) can be viewed as determining X w (t) in the limit 
t — > —00, provided the last term vanishes as t — > — 00 (this last property, however, may 
not hold for all solutions). The form (17) of the solution for the reservoir gives, in place 
of (11), the g-oscillator equation 

-t 



q + u q 



duo 







dt'qit') 



a 2 (co) 



00 



sin [00 (t — t')} e 



-o+\t-t'\ 



dcoa(u) < - [Ar(lo) + iBr(u)] exp(-iut) + c.c. > = 0. 



In the frequency domain, (17) is 

X u (u') = TT [A R (u) + iB R (u)] 6(0O - fa/) + 7T [A R (u) 



a{u)q{u') ma{u)q{uj') 

+ p — w + wr. i 6 ^ 



uo 



2co 



UJ 



iB R (uo)]5(u + u') 
- 8{uo + a/)] , 



(19) 



which inserted into (14) gives the following equation for q(to') 

i7ra 2 (|fa/|)" 



/2 



00 



doo- 



a 2 (u) 



,/2 



+ 



q(uo') 



uj' - bo'" 2uo' 

= -Tradw'l) [A R (\co'\) +isg a (co')B R (\co'\)} . (20) 

This is the frequency- domain version of ( |I8[ ). The solution q(co') of (20) gives the 
reservoir solution X u (t) by substitution of q(t) into (17). This second formulation of the 
dynamical equations presents a much easier path to solutions, as we avoid the integral 
equation (16). Moreover, as every solution can in principle be obtained by solving either 
(16) or (20), it may seem that (16) should be avoided entirely. But the ability in (16) 
to impose directly conditions on the reservoir at a finite time allows the discovery of 
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interesting particular solutions that in practice would not be found from (20). This will 
be clearly demonstrated in sections [3] and |4| 

Before proceeding to quantize the system for a general coupling function a(u>), 
we show classically that a particular coupling function gives damping proportional 
to velocity. This coupling function will therefore allow an exact quantization of the 
dynamics (JI]). 

3. Coupling for damping proportional to velocity I 



We consider the coupling function 

27 



a(ui) = uiquj 



-,1/2 



7 >0, 



(21) 



7r(w 2 + 7 2 ) 

which increases monotonically from at oo = to an asymptotic value of ujq^J as 



to — > oo. The function (|2l|) satisfies the integral relation 
P 



d^ 



u 2 + 7 2 



(22) 



and for this choice of coupling the integral equation (16) reduces to 
q(0 



d£ 



7TQU 

7T 
U 



[u 2 + 7 2 - ul] q\ 



2vr 

7 



(co 2 + 7 2 ) 



1/2 



[sgn(u)A (\u\)+iB (\tu\)}. 



(23) 



The general solution of (]23|) will contain the solution of the homogeneous equation 

7TU 



g(0 

e-o; 



7^o 



[w 2 + 7 2 -w 2 ] 



(24) 



We first solve the homogeneous equation (24) and then solve the inhomogeneous 



equation ( |23| ) for arbitrary v4 (u;) and B (u) by solving an associated integral equation 

the general solution of the homogeneous 



for a type of Green function. 

Following the standard procedure 



equation (24) is found to be 
q(u) 



(u 2 — Uq) 2 + 7 2 w 2 ' 



(25) 



where a is an arbitrary constant. Transforming (25) to the time domain, we obtain 



three different behaviours depending on the size of the constant 7 > 0: 



?(*) 



6exp 



7l*l 



COit\ 7 / U)i\t\ 

cos| T'] + ^ Bm IT" 



6exp(-wo|*|) (l + w |t|) , 



cexp 



7l*l 



cxp 



Ti 1*1 



7-7i 
7 + 7i 



exp 



Ti 1*1 



for 7 < 2w , (26 a) 
for 7 = 2u , (266) 
for 7 > 2u , (26c) 
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where b and c are arbitrary constants, and uj\ and 71 are defined by 



UJ L 



Acj 2 



r 



7i 



r 



Aoj 2 . 



(27) 



Note that u\ in (26a) is real, as is 71 in (26c). The presence of the absolute value of t 
in the solution (26a)-(26c) may at a glance appear to give a discontinuity in q(t) and a 
delta function in q(t), but in fact q(t) and q(t) are continuous for all t. 



The g-oscillator dynamics (26o)-(26c) is exactly the solution of 

0, 



g + 7g + u q 



0. 



t > 0, 
* < 0, 



(28) 



Q-1Q + u q 

with the condition g(0) = 0. We have thus obtained Ohmic damping and amplification 
of the g-oscillator in a closed system that can be canonically quantized. Recall that 
the time t = 0, which divides the damping and amplification epochs in (28), was an 
arbitrary choice. If we wish to impose arbitrary values of q and q at some time ti, with 
Ohmic damping for all later times, we need only replace t = in our derivations by a 
time to < ti, such that to an d the arbitrary constant in q(t) give the required values of 
q and q at ti. 

The dynamics (26a)-(26c) is the complete solution for the g-oscillator when (23) is 
valid and A (oj) = B (u) = 0. This last condition sets the displacements and velocities 
of all the reservoir oscillators equal to zero at t = (see ([6])). The solution for the 
reservoir is obtained from (|6j) by substituting (26a)-( p6~c ) and A (u) = B (u) = 0; the 
result is 



7 < 2a;o 
X u (t) = 

+ exp 
7 = 2w 
X u (t) = 

7 > 2a;o : 
X u (t) 

x 



buiuiQ 



2 7 



uj 2 — Uq) 2 + ^u 2 tt(u 2 + 7 2 ) 



1/2 



{ujq — 7 2 — u 2 ) cos(u;t) + 



(u 2 + 7 2 



Uq) COS 



1/2 



LO\t 



7 / 2 

H (uj 



7 



3lo ) sin 



7^o 

UJ 

Uli\t\ 



sin(a;|t|) 

(29a) 



[tt{u 2 + Au 2 ) 



exp(— oj \t\) [uj 2 + 3^0 + (w 2 + cj^^oltl] 



(w J + 3wo) cos(a;t) + 2— sin(a;|t|) 



(296) 



(uj 2 - uj 2 ) 2 + ^uj 2 L8vr(u; 2 + 7 2 ) 



7 



1/2 



571 



7^o 



(cJq — 7 — w ) cos(a;t) H sin(a;|t|) 



7 + 7i 

(4w 2 + (7 + 7l ) 2 ) exp 



7iN 



7-7i 



+ exp 



tI*I 



(4w 2 + (7 - 7l ) 2 ) exp 



li\t\ 



(29c) 



7 + 7i 

The under-damped case (26a) and (29a) is plotted in Figure [l] for a choice of 
parameters. Viewed from the infinite past to the infinite future, the g-oscillator is 





Figure 1. The solution (26a) and ( |29a[ ) with ujo = 3, 7 = 1 and 6 = 1. This solution 
is for coupling (21 ) with X u (0) = X u (0) = 0. The g-oscillator is damped into the past 



and future from t = 0, with damping proportional to velocity as shown in (28). As 
\t\ — > 00, most of the g-oscillator energy goes into reservoir oscillators with frequencies 
around ojq(= 3). 



initially (t — > —00) at rest with all of the energy in the reservoir; the g-oscillator is then 
amplified by the reservoir until at t = it has extracted all of the reservoir energy; 
then the energy is returned to the reservoir as the g-oscillator is damped. Alternatively, 
viewed as a consequence of the imposition of a t = condition, the g-oscillator is given 
a displacement, with g(0) = and X w (0) = X w (0) = 0, and is then damped into the 
past and future, transferring its energy to the reservoir. Note in Figure [T] that as the 
reservoir is amplified increasing from 0) and removes energy from the g-oscillator, 
the reservoir oscillators that attain the largest amplitudes as \t\ — > 00 are those with 
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frequencies close to uo, the free oscillation frequency of the g-oscillator. Compared 
to Figure [TJ the over-damped case 7 > 2u has the following qualitative differences: 
the g-oscillator is exponentially damped into the past and future from t = without 
oscillating, and the amplitudes of the reservoir oscillators for \t\ — > 00 decrease as u> 
increases from (so there is no peak in the reservoir amplitudes at u = u>q). 

The homogeneous version of equation (11) (i.e. with X u (0) = X w (0) = 0) has in 



fact a slightly more general solution q(t) than (26a)-(26c) for coupling (21). We did not 



obtain this more general solution in the frequency domain by solving the homogeneous 



integral equation (24), because its Fourier transform q(u) is not well enough behaved to 



make the derivation of (24) valid. For coupling (21), the time-domain equation (11) is 



q + u%q — sgn(t)7a;o / dt' g(t') exp(— j\t — t' 
Jo 



!) = /(*), 



/(*) 



du u u 



7 



_2tt(o; 2 + 7 2 )_ 



1/2 



(30) 



{[A (cj) + iB (u)] exp(— iut) + c.c.} 



which is (23) in the time domain. The homogeneous version of (30) (/(£) = 0) has a 



solution (26a)-(26c), but the more general solution is 
q(t) = Ub 2 - b x ) 



+ exp ( - 

+ exp 



2 

bi) 



b\ cos 



62 cos 



2 ) 



hi 2 + (62 - 



0J\t 

~2~ 



7^1 



b2l 2 - (p2 - frlVo 

7^i 



sm 



Uit 

~2~ 



sm 



CUit 



t > 0, (31a) 



t < 0, (316) 



where 61 and 62 are arbitrary constants and u± is again given by (27). For simplicity we 



have written (31a)-(31&) in the form appropriate for the underdamped case 7 < 2u ; 

2uq can be obtained by 



this solution is also valid for 7 > 2uo and the solution for 7 
taking the limit 00 1 — > 0. 



The solution (26a)-(26c) corresponds to the choice b 2 



in (31a)-(316). The more general solution (31a)-(316) differs from (26a)— (26c) by not 
having q(0) = and by having a constant displacement of the g-oscillator as \t\ — > 00. 



Note that the displacement q(t — > 00) in (31a) is minus the displacement q(t — > —00) in 
pill ). In the solution (glj-glj} the reservoir extracts all the kinetic energy from the 
g-oscillator as |t| — > 00 but does not bring it to its uncoupled equilibrium displacement 
q = 0. By inserting (31a)-(316) in ^ with A (co) = B (co) = 0, the more general 



version of the solution (29a)-(29c) for X u (t) is obtained, but we refrain from writing 
the result. 



To complete the solution for the dynamics with coupling (21), we must solve the 



equations of motion for non-zero X w (0) and X^O). We can write the solution in closed 



form for any Aq(u) and Bq(u) by finding a Green function for (30). If we find a solution 
G(t,t ) of 



G(t, t ) + oo 2 G(t, t ) - sgn(t) 7 ^ f dt' G(t', t ) exp 

Jo 



(- 1 \t-t'\) = S(t-t ), 



(32) 
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then the general solution of (30) is 



?(*) 



dt G(t,t )f(t ) 



(33) 



plus the solution (31a)-(316) of the homogeneous version of (30) (with f(t) = 0). The 
solution for X u {t) is then obtained from ([6]). A Green function G(t, t ) can be found by 



Fourier transforming (32) in t and solving the resulting integral equation for G(ui : to). 



As (23) is the Fourier transform in t of (30), we see from (32) that the integral equation 
for G(u, to) is 



TXUJ 
1^1 



[ w 2 + 7 2_ w 2] G(W;to) 



ti(u 2 + 7 2 ) 



exp(iut ) 



(34) 



Z-u 7^5 L ' 1Uq{uj - le) 

where a pole at u = in the last term has been moved off the real line by inserting 



e / 0. The two choices of sign of e in (34) give two Green functions G(t,to) in the 



time domain that differ by a solution of the homogeneous version of (32) (with on 



the right-hand side); this homogeneous solution is of course given by (31a)-(316) with 
a particular b% and b<i- The Green functions have rather lengthy expressions and the 



results for G(u,t ) and G(t,to) are given in Appendix A for the choice e > 



Although the classical dynamics is now completely solved for the coupling (21 ), the 



diagonalization of the Hamiltonian (in both the classical and quantum cases) will make 



use of an alternative formulation based on (17). 



4. Coupling for damping proportional to velocity II 



The forms (17) and (18) of the dynamical equations can be obtained by pushing back 



the (arbitrary) integration limit t = in ^ and (11) to t — > — oo. The part of X u (t) 
in (17) that is the solution of the homogeneous version of ^ (with q(i) = 0) describes 



the state of the reservoir at t — > — oo if the last term in (17) vanishes as t — > —oo. 



For the coupling (21), which satisfies (22), the g-oscillator equation (20) reduces to 

CO 



—(u +i'yu-u )q(u) = g{u), 
oj + 17 



-UUJq 



27T7 

uj 2 + 7 2 



1/2 



[sgn(w)A fl (M) + ifl fl (M 



(35) 



In the time domain this is (18) with coupling (21): 
q + u 2 q - 7^ 2 [ dt' q(t') exp[- 7 (t - *')] = 9®, 



9(t) 



7 



ll/2 



2n(u 2 + 7 2 



{[A R (u) + iB R (u)] exp(-iwt) + c.c.} 



(36) 



(37) 



Because of the factor of oj on the left-hand side of (35), the homogeneous equation 
(g(oj) = 0) has the solution 

q(u) = 2na5{uj) q(t) = a, (38) 
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where a is an arbitrary constant. The corresponding solution for the reservoir is obtained 
from (17) with q(t) = a and Ar(u) = Br(oj) = 0; the reservoir is also independent of 



time and the full solution is 



?(*) 



x u {t) 



27 



1/2 



(39) 



u itt(u 2 + 7 2 ) 

This "zero-mode" solution has zero energy, as can be verified using ([3]). It is interesting 
to check when a zero mode occurs in the general dynamical equations (El) and ([5]). 
Inserting q(t) = a into Q and ^ one finds that for o / there exists a consistent 
solution 

a(co) 



x u (t) 



if and only if the coupling function satisfies 

a 2 (co) 



0J n 



(40) 



(41) 



The coupling (21 ) is one example of a function that satisfies (41), and this is why there 



is a zero mode (39). Zero modes are always a possibility for oscillators coupled as in ([2]); 
even for just two oscillators coupled in this fashion there exists a zero mode for one value 
of the square of the coupling constant. In the case of the system ^ there is an infinite 
class of coupling functions that give a zero mode, and an infinite class that do not. It 
is interesting that the coupling function that gives Ohmic damping of the g-oscillator 
also allows a zero mode, but the wider significance of this fact is not immediately clear. 



Note that the zero mode (39) is not related to the constant displacement as t — > ±oo 
of the g-oscillator in the solution (31a)-(316); that solution has X w (0) = X w (0) = 



and so has no relation to (39). The formulation of the dynamical equations used in this 



section has shown that the zero mode (39) is completely determined by the condition 
Ar(u) = Br(oj) = on the reservoir, and we see from (39) that this is not equivalent 
to the condition X w (t)| t _ ) ._ 00 = 0. 



The general solution of (35) is the solution (38) of the homogeneous equation plus 
the solution 



q[u) = 
G R (u) 



(w + i0+)(o; 2 + i7^ 



(42) 
(43) 



In (43) we have defined a Green function and given a prescription for dealing with its 
pole at uj — 0. We have moved the u = pole into the lower-half complex w-plane, where 
the other two poles of Gr(uj) lie, and have therefore chosen the retarded Green function 
(analytic in the upper-half u- plane). It is easy to see that a different prescription for 
dealing with the u = pole will give a Green function that differs from (43 ) by a solution 



(38) of the homogeneous equation. The Green function (43) in the time domain is 

2 



G R (t) 



m- 2 



-0+t 



exp 



jt 
' 2 



cos 



+ 



7 



2u$ 



sm 



7^i 



UJ\t 



(44) 
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with oj\ given by (27). The result (44) is valid for all values of 7 if the case 7 = 2cuq 



is understood as the limit u)\ — > 0. As is expected from (36), the Green function (44) 
satisfies 

(45) 



G R + <4G R - 7^ 2 / dt' G R (t') exp[- 7 (t - 0] = 6{t). 

J —00 

In the time domain the general solution for the g-oscillator is 

/>oo 

q(t) = a+ / dt'G R {t-t')g(t>), 



(46) 



as can be confirmed by showing that the second term satisfies (36) because of (45). The 



solution for X u (t) is given by (17) 



By removing the step-function factor 6{t) in the retarded Green function (44), we 



obtain a solution of the homogeneous version of (45) (with zero on the right-hand side) 



This is of course also a solution of the homogeneous version of (36), so there seems 
to be a more general solution for q(t) with A R (u>) = B R (u) = than the zero-mode 



solution (38). In fact the general solution of the homogeneous version of (36) is given 



by ( 31a[ ), with the expression valid for all t rather than t > 0. (The Green function 



(44) with 6 {t) removed corresponds to (31a) with b\ = — 7/w 2 and b 2 = j/uq.) But the 



expression in (31a), taken to be valid for all t, is not a solution for the dynamical system 



for the simple reason that the corresponding reservoir solution, given by (17), diverges. 



This behaviour of the dynamics can be described with reference to Fig. [TJ which shows 



a particular case of the solution (31a)-(31&) and the corresponding reservoir solution. 
The time t = in Fig. [I] is arbitrary and the solution holds with the t = peak in q(t), 
and zero in X u (t), moved to any other finite value of t. If this peak in q(t) is pushed 
back to t — > —00, however, the amplitude of the peak diverges and this causes the entire 
solution X u (t) to diverge. The general solution for the dynamics with the condition 



B r (lu) = is therefore the zero mode (39) 



The solution shown in Fig. [T] must be a particular case of the general solution as 
formulated in this section: it must correspond to some choice of A R (cu) and B R {u). As 
the solution in Fig. [ljhas q(t) — > as \t\ — > 00, the corresponding functions A R (co) and 
B R (u) turn out to specify the state of the reservoir at t — > —00 (see (17)). By taking 



the limit t — > —00 in (29a) we see that X u (t) takes the form of the first two terms in 



(17), with 



A R (u) 



B R (u) 



-b 



2 7 



1 1/2 ww (w 2 + 7 2 - ^o) 



ir(u 2 + 7 2 ) 
27 

ii(u 2 + 7 2 ) 



-,1/2 



(u 2 — Uq) 2 + ^ 2 0J 2 
(OO 2 — U)q) 2 + 7 2 Co> 2 



(47) 



(4* 



One can verify that with A R (u) and B R (cu) given by (47) and (48), and with a = 



equations (46), (44), (37) and (17) reproduce the solution (26a)-(26c) and (29a)-(29c) 



Needless to say, it is by far from obvious that the specification (47) and (48) of the 
reservoir at t — > —00 implies the interesting behaviour depicted in Fig. [l] (modulo a 



zero-mode solution (39)). 
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We have formulated the general solution for the dynamics with coupling (21 ) in two 



ways. Every particular solution is completely specified by either the functions Aq(u) 
and Bq(u) and constants 61 and 62 in section |3j or by the functions A r (oj) and B R (tu) 
and constant a in this section. The relationship between the pairs {Aq(u), Bq{oj)} and 



{A r (uj), B r (uj)} for the same solution {q(t), X u (t)} is the one that transforms (13) into 



(19) and (16) into (20); it is easy to show that this relationship is 



a n n a n in "(M) f T r / \ 1 2u; f°° , Re[g(011 
A {\u\) = A R (\u\) - lm[q(u)] + — P / d£ } 



2u 



71 



S 



Bo(\u\) 



B R (\u\) + sgn( W )^M J Re[q{uj)] _ -p 
2a; 7r 



d ^Im[g(e)] l 
^ 2 — 



■ 



(49) 
(50) 



The results (47) and (48) for the particular solution ( |26a[ )-( p6c ) and (29a)-(29c) can 
also be derived from (49) and (50) by inserting A (co) = B (cj) = 0, (21) and (25). 



5. Quantization and diagonalization of the Hamiltonian 

The canonical momenta for the Lagrangian (|2j) are 
Tl q (t) = q(t), U Xui (t)=X u (t). 



(51) 



H 



2^ + 2 W ° g + 2 



dui a(u) 



qX^ + X u q 



(54) 



We quantize the system in the Heisenberg picture by imposing the equal-time canonical 
commutation relations 

[<?(*), n,(*)l [X UJ (t),fl Xbj ,(t)} = ih6(u-co'), (52) 

[X w (t),X w ,(t)] = 0, [fl Xu (t),fl Xujl (t)} = 0, [q(t),X u (t)] = 0, [g(t),fl x Jt)] = 0.(53) 

The commutation relations for the reservoir are similar to that of a field theory because 
of the continuum of frequencies u. The Hamiltonian is 

where a Hermitian combination of the operators has been taken in the last term; 
classically, (54) gives the total energy 

The solution of Hamilton's equations for the canonical operators can be immediately 
written down using the classical results in sections [3] and [4} But in order to describe 
general quantum states we need some basis in the Hilbert space, and a solution for 
the canonical operators does not by itself reveal such a basis. Diagonalization of the 
Hamiltonian solves Hamilton's equations for the canonical operators in terms of the 
energy eigenstates of the system and thus solves the dynamics of general quantum states. 
The diagonalization process can also be performed on a purely classical level, where the 
classical normal modes of the system correspond to the quantum energy eigenstates. 
What is required is a new set of dynamical variables for which the system reduces to 
a set of uncoupled harmonic oscillators; each of these free oscillations (normal modes) 
has a conserved energy and gives an energy eigenstate in the quantum theory. 

For a discrete set of oscillators with coupling terms of the same form as in ([2]), the 
normal-mode displacements are linear combinations of the displacements of the coupled 



Quantum dynamics of the damped harmonic oscillator 



16 



oscillators [19] . The continuum reservoir in ([2]), however, has the remarkable effect that 
the normal modes are given by a canonical transformation of the original variables, in 
which displacements and canonical momenta are mixed. It is well known that in the 
discrete case the frequencies of the normal modes become complex when the coupling 
of the form ^ is too large [T5] . We shall also find that there is a restriction on the 
coupling function a(oo) in order for the Hamiltonian (54) to be diagonalizable with real- 
frequency eigenmodes. This restriction on a(oo) is met by the particular coupling (21). 
In contrast, Huttner and Barnett [21] showed, as part of their model of a dielectric, 
that the system Q with in the coupling term replaced by has real-frequency 
eigenmodes for essentially any coupling function a(uo). 



We attempt to show that the Hamiltonian (54) can be written 
1 



H 



doo 



TT 2 



+ oo 2 $ 2 



doo hoo 



t)C{oo, t) + C{oo, t)C\oo, t) 



(55) 
(56) 



where $ w (i), oo G [0, oo) are displacement operators for a set of uncoupled harmonic 
oscillators, the normal modes or energy eigenmodes of the system. The normal-mode 
creation and annihilation operators C\oo, t) and C(uo, t) are given by the usual harmonic- 
oscillator expressions 



C(oo,t) 




00 



(t) 



C\oo,t) 




oo 

2h 



$ w (*)--n* w (t) 



00 



and they obey the commutation relations 

[C(oo,t),C\oo',t)] =5(oo-u'), 



[C(oo,t),C(oo',t)] = 0. 



From (56) and (58) we obtain 

hooC(oo, t), 



(57) 



(5f 



(59) 



[C(oo,t),H] 

so that the eigenmode creation and annihilation operators have a stationary time- 
dependence given by 

C{oo, t) = exp(-ioot)C{oo). (60) 



The diagonalized form (55) and (56) of the Hamiltonian (54) will be achieved with a 
restriction on the coupling function a{oo). The following derivation can be performed 
classically if commutators are replaced with Poisson brackets. 

If the diagonalization is possible then there exists a transformation from the original 
dynamical degrees of freedom in (54) to eigenmode degrees of freedom for which the 
Hamiltonian takes the form (55). In the quantum theory it is convenient to use the 
creation and annihilation operators (57) as the eigenmode degrees of freedom, rather 
than the canonical operators l?^ and 11$^ (the quantities (57) can also be used in the 
classical theory). We thus seek a transformation 



doo 



f q (u)C(u,t) +h.c. 



n g (t) 



doo 



f n (u)C(uo,t) +h.c. , (61) 
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da/ 



du/ 



fx(u,u')C(u',t) + h.c. 
f Ux (u,co')C(co',t) + h.c. 
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(62) 
(63) 



for which (54) becomes (56). The unknown coefficients f q (cu), etc., in (61)-(63) can be 
written as commutators by utilizing (58): 

f q (u) = [g(t),CV,t)], /n» = [n ? (*),C't(a;,t)], (64) 
f x (u,u') = [X u {t),&{u/,t)\, /n x (u;,u/) = [n^(0,CV,0]- (65) 

The eigenmodes variables must be expressible in terms of the original variables, i.e. 
(61)-(63) must be invertible, and (64), (|65j), (52) and (53) imply 



(66) 



C{uj,t) = --\r Uq (uj)q(t)-rMTl q (t) 



+ 



du' 



r nx {j^)x ul {t)-r x {u',u)u Xut {t) 



It is clear from (61)-(63) and (60) that the /-coefficients are closely related to 
Fourier transforms of the original canonical operators. This suggests that the f- 
coefficients must satisfy the original dynamical equations written in the frequency 
domain. We derive these equations for the /-coefficients as follows. By inserting (66) 
and (54) in (59), and using the canonical commutation relations (52) and (53), we obtain 

hw6(u>,t) = f^ q (u)tl q (t)+u 2 r q (u)q(t) + y°°du/ {f^ x (co',u)U x Jt) 

+u' 2 f x (cu',u)X UJ ,(t) - a(u/) f/!(w)t(i) + f x (u/,u)q(t)} } . (67) 



Comparing coefficients of the canonical operators in (66) and (67) we find 

POO 

fn q (u) = iw/n,(w) =Wo/,(w) - / du'a(u')f x (u', w), 



fn x (u',u) = -hjfx(w',w), iufn x (w',u) = u' 2 f x (u', w) -a(u')f q (u), 
which give 

POO 

u 2 f q (u) = ulf q (uj) - / du)'a(u)')f x (u)',U)), 



u 2 f x (u',u) = u' 2 f x (u',u) - a(u')f q (u). 



(68) 
(69) 

(70) 
(71) 

These two equations for f q (cj) and f x (co', u>) are indeed identical to the frequency-domain 
equations (14) and (15) for q{u) and X u i(cu). We can therefore write the general solution 
of (70) and (71) using the results (19) and (20), in which a retarded Green function was 
used to solve the reservoir equation. As only positive frequency arguments are used in 
the /-coefficients, we can drop delta functions containing sums of frequencies and so 
(19) gives the following general solution for f x (u',u): 



f x (u',u) = h x (u)5(u-u') + P- 



2u' 



a(cu') 



(72) 
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where h x (u) is an arbitrary complex function. The corresponding general solution for 



f q {oS) is found from (70) and (72), which yields (20) with positive frequency arguments: 



u, ! + P 



s C2 



+ 



ma (uj) 



2cu 



fq{u) = -a(u)h x (u). 



(74) 



The solution of (74) can be written 



f q (ui) = h q {u) + G(ui)a(u)h x (u) 



G(u) 



UJ 



« 2 (0 



+ 



i7ra I a; 



(75) 



where we have defined a Green function G(u), and h q 



.u>. 



2uj 

is the solution of 



UJ' 



+ 



17TOJ {UJ) 



2w 



= 0. 



(76) 



The existence of a non-zero solution of ( 76 ) requires at least one real root of the function 



in brackets; the resulting h q (u) will contain a delta-function factor and an example is the 



zero mode (38) for the coupling (21). The /-coefficients are now completely determined 
by the two functions h x {u) and h q (uj); it remains to show that values for these two 
functions exist that give the required diagonalization. As follows from section [2j there 
are alternative ways of writing the general solution for the /-coefficients, but it turns out 



that the form (72)-(76), based on the retarded solution of the reservoir equation, is very 



convenient for solving the diagonalization problem. An equally convenient approach, 
which is in fact adopted by Huttner and Barnett [21J, is based on the advanced solution 
of the reservoir equation. 

A value for the function hx{w) is obtained by demanding that the commutation 



relations (58) hold for the operator (66) and its hermitian conjugate. The first of (58) 



after elimination of /n„(^) and fjj x (u',u) using (68) and (69), yields 



{uj + u')r(uj)f q {u') + {uj + u') / dco" r x (u",u)f x (u",u') = h5(u - to'). 



(77) 



The key part of the diagonalization is the insertion of (73) into (77). A product of 



terms containing the infinitesimal number + is encountered, whose value depends 



quadratically on + . This means that the form (73) of f x (uj',uj) must be used in (77) 



rather than the form (72). Any diagonalization involving the continuum reservoir will 



encounter such products, and their treatment is described in detail in [22J, where the 



Hamiltonian of macroscopic QED is diagonalized. The result of substituting (73) into 



(77) and applying (74) is 

2uh* x (oj)h x (u) 
which has a simple solution 



h. 



h x (ui) 



-) 



1/2 



(78) 



(79) 



Other solutions of (78) differ from (79) by a phase factor. This freedom in the choice of 



the diagonalizing transformation is clear from the outset, as the diagonalizing operators 
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(57) are only defined by (56) up to a phase factor. A similar analysis shows that the 
second commutation relation in (58) is identically satisfied by (66) because of (73) and 



We must also check that the commutation relations (52) and (53) of the original 
canonical operators are satisfied by the representations (61)-(63). This is a consistency 
check on the /-coefficients, which have already been largely determined by ( 78 ) . It is here 
that we find a restriction on the coupling function a(u). The full set of coupling functions 
that satisfies the restriction is determined by rather complicated integral relations, but 
for a subclass of coupling functions that includes (21) the restriction is much simpler. 
Consider the set of functions a(u) for which 



a (to) is an even function of to 



and 



to 



> except possibly at to — 0. (80) 



For coupling functions satisfying (80) the diagonalization can be achieved if and only if 



either lo > 



or to. 



o 



and uj n + 



d£ 



« 2 (0 

£ 2 + K 2 



(81) 



~ K n , n < 1, as k — y 0. 



The proof of this is given in Appendix B , which furthermore shows that the diagonalizing 
transformation requires the choice 



h q (to) 







12) 



as the solution of (|76j). It is also shown in Appendix B that for the class of coupling 
functions (80), there is only one possible real root of the function in brackets in (76), 
at to = 0. Moreover this to = root occurs if and only if ujq = / °° d£a 2 (£)/£ 2 , which 
we recognize as the condition (41) for the existence of a zero mode. Thus, the second 
condition in (81) is a restriction on coupling functions that have a zero mode. 

Assuming the conditions (80) and (81) are met, a set of /-coefficients for the 
diagonalizing transformation is specified by (73), (75), (79) and (82). The resulting 
relations (61)-(63) for the canonical operators in terms of the eigenmode creation and 
annihilation operators are most simply written in the frequency domain. With the 
Fourier definition 

1 



2tt 



(to) exp(— kot) + h.c] , 

for q(u), etc., the frequency-domain canonical operators are 

h\ 1/2 



2tt 



q(to) = 27r 



2oo 



- I a(co)G(io)C(io 
a(cu) 



2uo J 

5{u-u')C{uj') + 



-U q (u), 



1 



2u \ uo — to' — i0 + oj + uj 1 



q(u f ) 



UJ' 



(83) 

(84) 

(85) 
(86) 
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where the Green function G(oS) is defined in (75). The Green function can be written 

-1 



G{u) 



UJ A — Un 



[i-xH] 



(87) 



in terms of a dimensionless quantity x( w ) that gives the dependence on the coupling 
function a(u): 



« 2 (0 



e-^ + ^r- m 

is an even function of u then x( w ) exhibits a Kramers- 
Kronig relation between its real and imaginary parts and is therefore analytic in the 
upper-half complex cu-plane. We can view as an effective "susceptibility" for the 
damped g-oscillator; remarkably, for couplings such that a 2 (oj) is even, a single damped 
oscillator exhibits dispersion and dissipation connected by Kramers- Kronig relations. 
There are in fact many similarities between the quantum damped harmonic oscillator 
and macroscopic QED (22], where the incorporation of Kramers-Kronig relations for 
electromagnetic susceptibilities is an essential part of the theory. 



As noted in Appendix B if a (u 



In sections 3 and 4 we solved the classical dynamics for the coupling function (21). 



We note that this coupling function is of the class (80). In addition, the function (21) 
has the properties 

■« 2 (0 _ 2 7 c 2 



Z i z 
CJn + K - 



d£ 



£ 2 + K 2 



+ r 



CO + K 



o- 



1^1 

K + 7 



7 



(89) 
(90) 



The second of the restrictions (81) is thus satisfied, and therefore the Hamiltonian is 



diagonalizable in this case. The effective susceptibility for the coupling (21) is, from 



(88), 



x(w) 



7 



7 



The Green function (87) is then 
G R (u) = 



2\ ' 



(91) 



(92) 



(to + i0+)(o; 2 + i'jtu — co 0) 
where we have moved a pole at u = into the lower-half complex plane; this gives the 
retarded Green function Gr(u), analytic in the upper-half plane. The retarded Green 



function (92) is identical to (43), which featured in the classical solution (as discussed 
above, the /-coefficients are a particular solution of the classical dynamical equations). 
The choice of the retarded Green function is not required for the diagonalization; other 
choices will give /-coefficients that differ by a zero-mode solution, described in section |4j 



6. Coherent states 



Having diagonalizaed the quantum (and classical) Hamiltonian for coupling functions 
that include (21), we can now construct the quantum state of the system that is closest 
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to the interesting classical solution derived in section [3] and plotted in Fig. [Tj We expect 
the quantum states that are closest to classical solutions to be coherent states, and one 



benefit of the diagonalization results (84) and (85) is that they allow us to construct 



coherent states of the eigenmodes of the system. Classically, what we are about to do is 
to relate the particular solution in section |3]to the classical normal modes of the system. 

The eigenmode annihilation operator C(u) defines continuous- mode coherent 
states [M] by 



C{u)\C{u))=C{u)\C{u)), 



(93) 



where C(cu) is an arbitrary complex function. The classical solution for coupling (21) 



gives the relations (35), (42) and (43) for q(uj), while the quantum solution for q(oo) with 



coupling (21) is given by (84) with Green function (92). Comparing these classical and 



quantum results we see that every classical solution given by the real functions A R (ui) 
and B R {u) has a corresponding quantum coherent state with complex amplitude 

C(u) = -[^) 1/2 [A R (uj)+iB R (u)]. (94) 

The expectation values of q(t) and X w (t) in the coherent states are equal to the 
corresponding classical solutions q(t) and X u (t). It follows from the results at the end 



of section H that the coherent state corresponding to the classical solution (26a)-(26c) 



and (|29a|)-([29c|) has amplitude Q with A R (u) and B R (u) given by (|47[ and Q. 
The expectation value of the normal-mode displacement operator $ w (t) in a 



coherent state is, from (57), (93) and (94) 



(C(u)\^(t)\C(u)) 



( ft \ 1/2 

( — j \C{uS) exp(— iutf) + c.c.] 



(95) 

= — A r (oj) cos(a>t) — B R {u) sin(u;t). (96) 

This is just the general solution for the classical normal modes ^ u (t) of the system. 
There are many other results that can be derived for the coherent states, but in this 
paper the quantum state we will consider in detail is the thermal state. 



7. Thermal equilibrium 

Another benefit of diagonalizing the Hamiltonian is that it allows a straightforward 
calculation of the thermodynamic quantities for the g-oscillator. In this section we 
treat the case of thermal equilibrium for all couplings for which the diagonalization in 



section 5 is valid; the results for the particular coupling (21), corresponding to damping 



proportional to velocity, will also be described. Given the close relationship between the 
damped harmonic oscillator and macroscopic QED, noted in the last section, it is not 
surprising that the following treatment of the thermal state of the damped oscillator 
has similarities to the Casimir effect, which is simply macroscopic QED in thermal 
equilibrium [23J. 
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In thermal equilibrium the excitation level of each eigenmode oscillator of the 
system will be given by the Planck distribution, so the thermal mixed state is defined 
by 

(C\co)C{co')) = M{co)5{co - co') = {C{co)6\co')) - 5(co - co'), (97) 

(98) 
(99) 



N{u) = exp 

\k B T 
(C{co)C{co')) = 



The continuum of eigenmodes necessitates the use of delta functions in (97). These 



delta functions lead to some awkwardness in the thermal calculations, but as in 
such difficulties can be negotiated by regarding such delta functions as limits to be 



strictly imposed only at the end of calculations. Using (97)-(99) it is straightforward to 



compute the thermal position and momentum correlation functions for the g-oscillator, 
as well as its thermal (including zero-point) energy and other thermodynamic quantities. 
Our procedure will be to derive the results for a general coupling function for which the 
diagonalization in section [5] is valid, and then to examine those results for the particular 



coupling (21 ). 



7.1. Thermal correlation functions for the q-oscillator 



It follows from (83), (84) and (99) that the thermal correlation function of q(t) can be 
written 

1 

4^2 _ 

(100) 



(mm) 



dco / dco' {exp [-i(ut - co't')} (q(w)tf(u/)) 
Jo 

+ exp [i(ut - co't')} (gt(w)g(w'))} , 
with a similar result for the other canonical operators of the system. The frequency- 



domain correlation functions (q* (u)q(w')) and (q(u)q* (u)')) are, from (84) and (97) 
*t/..w..m 27rW(u;) 



(<f(u,)g(u/)} 



CO 

M{oo) 



M{u) + 1 



G*{co)G{co')N{co)5{co - co') 
(q{u)q\J)). 



(101) 



Multiplying the Green function expression (87) by G*(co), we obtain 
{co 2 - co 2 [1 - x(u)]}G*{co)G(co) = -G*(co), 



and the imaginary part of this equation is found using (88): 



it a {co) 
2oo 



G*(uo)G(oo) =lmG(u) 



Inserting (103) in (101) we find 



(tf(co)q(co')) = 4ttM(w - co')Af(co)ImG(co) 



Afjco) 
jV(co) + 1 



(102) 



(103) 



(104) 
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which shows that the real part of the temporal correlation function ( 100 ) is 

h 



mm + mm) 



71 



dui cos [u(t — t')] coth 



2k B T 



(105) 



where 2j\f{oS) + 1 has been rewritten as a hyperbolic tangent. The momentum operator 
for the g-oscillator is given by (84) and is just the time derivative of q(t), so we obtain 
from ( 105[ ) the momentum correlation function 



fl q (t)fl q (t') + fl q (t')fl q (t) 



h 



7T 



dcu ijj 1 cos [u(t — t')] coth 



2k B T 



lmG(u). (106) 



We now consider the general results (105) and (106) in the case of the coupling 
function (21 ), corresponding to damping of the g-oscillator proportional to velocity. We 
first note a slight modification of the result (103) in the case of (21). Equation (102) 
does not include a possible pole prescription at u = for a zero mode, such as occurs in 
(92). Evaluating the left-hand side of (103) for the retarded Green function (92), and 
replacing + by rj, we find 



ixa [oj] 
2u 



7^o 



(107) 



{u 2 + r] 2 )[{u 2 - u 2 ) 2 + 7 2 w 2 ]' 

As 7] — > 0, the right-hand side of ( 107) has a pole at u> = and r] prescribes the treatment 
of this pole by a principal value; apart from the pole prescriptions, (107) is the result 
(103). With use of (107) the correlation function (105) gives the following expectation 
fit)) for coupling (21): 



value 



h 

2^ 



dui 



-fUU 2 



coth 



hw 



(108) 



(co 2 + r] 2 )[(oo 2 - io 2 ) 2 + -i 2 u 2 } ~"~ \2k B T / 

In (108) the integration has been extended over negative u (the integrand is an even 
function of u); this allows evaluation of the integral by closing the integration contour 
in the upper (or lower) half complex plane. It is clear that the integral in (108) does 
not converge when rj = 0, and evaluation by contour integration for 77 > shows that 
this divergence is given by a term k B T^i /(uqI]). The form of the divergent term shows 
that (108) is finite at T = 0, so that the ground state value of (q 2 (t)) is well defined. 
The divergence at T > can be attributed to the zero mode, as it is associated with the 
bj = pole in the Green function. As the zero mode has zero energy, it is perhaps not 
surprising that the zero-mode squared displacement diverges at non-zero temperature. 
Although it is tempting to drop the divergent zero-mode term ksT'j / '{u>q1]) from (108) 
when T > 0, this would be a violation of quantum mechanics because the canonical 
commutation relations of the system only hold for the full displacement operator q(t). 
There is thus an inherent pathology in the thermal state for the coupling (21 ), although 
the ground state is well-defined as are other quantum states such as the infinite class of 
coherent states described in section [6| Note that the divergent term k B T^ / {uj^ri) also 
vanishes if the limit 7 — > in the integration result is taken before the limit rj — > 0; in 



As 



this case the result for (q 2 (t)) reduces to the free-oscillator value 
discussed in section [5j coupling functions in the class (80) give rise to a zero-mode only if 



^coth 

2lu 



2k„T 
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they obey the special condition (41); for general coupling functions in the class (80) the 



correlation function (105) will be well behaved. The momentum-squared expectation 



value (Ilz(t)} for the coupling (21) is found from (106) and (107) to be 



h 

2^ 



do; 



If + 

where there is no longer any contribution of the u - 



coth 



(109) 

2k B T / l 

pole in the Green function. 
The integral in (109) is well defined and can be evaluated by contour methods but, 
given the pathology in the corresponding position quantity (q 2 (t)) at T > 0, the rather 
complicated details are omi tted here. We note however that when the the limit 7—7-0 
is taken in the final result, (109) gives the free-oscillator value |fta;ocoth f gtfr 



7.2. Reservoir contributions to thermal correlation functions 

The thermal correlation functions of the reservoir are of interest because of their 
contribution to the thermal energy of the g-oscillator. As discussed in section |7.3 



for this purpose we must calculate correlation functions containing the reservoir using 



(85) and (86), but only retain terms that contain q(u). This is the prescription used in 
macroscopic QED to obtain the Casimir stress-energy of the electromagnetic field in a 
medium, where the electric and magnetic fields play the role of the g-oscillator [23J. In 
fact the results in this subsection are largely a simpler version of the reservoir correlation 
functions in |23|. 



The expectation value of the Hamiltonian (54) in thermal equilibrium contains 



the term J °° do;((9 t X aj ) 2 + u 2 X 2 )/2. We compute this expectation value using (85), 
retaining only terms in which q(u) occurs. Essentially the same calculation is described 
in detail in [23] for macroscopic QED. We therefore refer the reader to [23J for the 
various steps in the derivation and here state the result: 

1 f°° i //o O X2 2,> 2 \ H 2 r r° , , / *U \dbxM] 



da; 



;no) 



The thermal expectation value of the Hamiltonian (54) also contains f °° do; a(u)(qX u ) 



This expectation value is also found using (85), but here q(u) occurs in all contributions 



to the expectation value, so no terms are dropped. There is again an essentially identical 
calculation in [23] for macroscopic QED, and the same derivation gives the result 



dua(w) (q(t)X u (t) 



Hu 2 



-Im 



7T 



du coth 



hu 
2k B T 



xMG(w). 



Ill) 



This quantity is real so the term containing (qX u + X u q)/2 in the Hamiltonian (54) gives 
the thermal expectation value on the right-hand side of ( jlll ) . The thermal (including 
zero-point) energy of the g-oscillator can now be obtained. 

7.3. Thermal energy of the q-oscillator 

The energy of the g-oscillator in any state is given by the expectation value of the 



Hamiltonian (54) with the reservoir X w traced out. Due to the coupling of q to X u 
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this tracing operation is in general not a trivial matter. The same issue arises in 
macroscopic QED, where the value of the electromagnetic energy-momentum tensor 
in a dispersive, dissipative medium requires the tracing out of the reservoir. In the 
case of thermal equilibrium the electromagnetic stress-energy was found in [23] to 
be given by a straightforward recipe: after the reservoir operators are expressed in 
terms of electromagnetic-field operators, all contributions containing electromagnetic- 
field operators are to be retained, the other contributions are to be dropped. We follow 
the same recipe here to trace out the reservoir from the total energy of the system in 
thermal equilibrium and thereby obtain the thermal energy of the g-oscillator. The 
reservoir operators are expressed in terms of the g-operator using ( 85 ) and ( 86 ) , and the 
thermal expectation value of the Hamiltonian (54) is evaluated with only contributions 
that contain q(u) included. The result is the thermal energy of the g-oscillator, which 
we denote by (H) q . The required expectation values involving the reservoir are given 
by (110) and (111), and the thermal expectation values of the first two terms in the 
Hamiltonian (54) are found from (105) and (106). This gives the following expression 
for (H) q : 



h_ 

47T J 



Tm 



du coth 



2k B T 



dco 



G(w). 



;ii2) 



Due to the parity properties of the susceptibility and the Green function, the real part 
of the integrand in (112) does not contribute and the integral is imaginary. The q- 
oscillator thermal energy (112) is the analogue for the damped oscillator of the Casimir 
energy density of electromagnetic fields in a medium [23]. A dispersive contribution to 
the energy from the derivative of the susceptibility is familiar in the electromagnetic 
case [23] , and we find a similar contribution from the damped-oscillator "susceptibility" 
in (112). As is familiar in Casimir theory and elsewhere, the integral in (112) can be 
converted to a sum over imaginary Matsubara frequencies. 

In the case of the coupling (21 ) with susceptibility (91 ), the real part of the integrand 



in ( 112 ) does not give a convergent integral so that it is not correct to take the Im outside 
the integration as has been done in the derivation of (112). For the coupling (21) we 
must keep the Im inside the integration and take the imaginary part of the integrand; 
the resulting expression for (H) q is 



h 

4tt 



du 



(to 2 + 7 2 ) [(u 2 - u 2 ) 2 + 7 2 u; 2 ] 



coth 



2k B T 



113) 



There is no zero-mode contribution in (113), as is expected from the fact that the 
zero mode has zero energy. The thermal energy (113) is finite and can be evaluated 
analytically by closing the integration contour in the upper (or lower) half complex 
plane. When the zero-damping limit 7 — > is taken, however, the result for the thermal 
energy is found to reduce to 



-hwo coth 



hw 
2k B T 



k B T, 

2 



;ii4) 
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i.e. in addition to the free-oscillator thermal energy there is an additional term —ksT/2. 
This additional term can be traced back to the dispersive contribution dx(w)/du in 
(112); although vanishes as 7 — > 0, the contribution from this dispersive term to 
the energy becomes — ksT/2 as 7 — > 0. The non-zero contribution of the dx(w) / dw-term 
when 7 — > relies on the assumption that the susceptibility (91) is valid up to infinite 
frequencies. Although it is mathematically convenient to employ functions like (91) over 
the entire frequency range, this is an unphysical assumption and therefore may produce 
some unphysical results. As will be discussed in section [8j application of the general 
results presented here to real physical systems would probably require a measurement 
of the susceptibility characterizing the particular damped oscillator; although these 
measured susceptibilities could be fitted to mathematical functions, non-zero values of 
these functions at very large frequencies beyond the measured range would be physically 
meaningless. We note that the difficulty with the zero-coupling limit of the thermal 
energy for damping (21) does not occur for the ground state T = 0, as we see from 
(114) that the correct free-oscillator limit huo/2 is obtained. These thermal-energy 
results support the comment in section 7.1 that the coupling (21) gives an inherent 
pathology in the non-zero temperature state. 



8. Conclusions 



The use of a discrete set of oscillators as a reservoir is almost universal in studying 
dissipation in quantum systems using exact quantization rules. Yet the power of a 
continuum reservoir in the study of dissipation was demonstrated twenty years ago by 
Huttner and Barnett [21], and it has recently been shown that a continuum reservoir 
allows an exact canonical quantization of the macroscopic Maxwell equations in arbitrary 
media obeying Kramers-Kronig relations [221 E31 EH E5]. The important results that 
have been obtained with a continuum reservoir would not have been achieved if a limiting 
procedure had been applied to results featuring a discrete reservoir, whereas the discrete 
reservoir almost never captures the desired physics without a subsequent, and delicate, 
limit to the continuum case |19j . In view of this, it is highly advisable to employ a 
continuum reservoir from the outset. We have considered a single damped quantum 
harmonic oscillator and shown that a continuum reservoir allows the exact treatment of 
damping proportional to velocity. In addition we have found that for general damping 
the quantum oscillator is remarkably similar to light in macroscopic media. For functions 
a(u) coupling the oscillator to the reservoir, with a 2 (u) an even function of u, the 
dynamics of the damped oscillator is governed by an effective susceptibility obeying 
Kramers-Kronig relations. 

The experimental investigation of macroscopic oscillators that exhibit quantum 
behaviour has made remarkable progress in recent years; oscillators have been cooled to 
their quantum ground states and even placed in a superposition of energy states [26l E7J 
[28| [29| [30] . A theoretical description of nano-oscillators and opto- mechanical systems 
can be approached in various ways, depending among other things on how much detail 
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of the microscopic structure of the oscillator is included. The approach to the damped 
oscillator taken in this paper naturally leads to a description of the oscillator by means 
of an effective susceptibility. As in macroscopic electromagnetism, the strength of this 
approach is that the complicated microscopic substructure is not included, rather the 
susceptibility is a quantity to be measured for real physical systems. We have found 
that diagonalization of the Hamiltonian requires restrictions on the coupling function 
that determines the susceptibility, but this may not be significant for real macroscopic 
oscillators. The condition (81 ) constrains the relationship between the coupling function 
a(u) and the free-oscillation frequency ojq. But most mechanical oscillators do not 
exhibit a "free" oscillation state; the oscillation degree of freedom is produced by the 
same material geometry that causes damping. This means that the parameter uq cannot 
be separated from the susceptibility, and they must together be fitted to a given system 
through experimental measurement. A choice of cuo for which Uq > / °° du a 2 (u)/u 2 may 
therefore not be restrictive. This also suggests that the there is freedom in the form of 
the coupling term between the g-oscillator and the reservoir. We chose a coupling of the 
form qX^, whereas Huttner and Barnett [21] use qX u , and the latter does not give the 



restriction (81) on the coupling function. But the Huttner-Barnett coupling leads to 
a renormalization of the free-oscillation frequency ujq |21| . which is why the restriction 
(81) is avoided. If ujq has no direct physical meaning one can effectively make the same 
redefinition with the coupling term chosen here, through the joint characterization of 
u and the susceptibility 

The considerations in this paper are of course only a first step towards describing 
physical quantum oscillators, and other ingredients would have to be included such as 
coupling to light in the case of opto-mechanical systems. There is also the crucial 
question of how the effective susceptibility is to be measured for real oscillators, 
or how physically realistic susceptibilities could be deduced from simple microscopic 
considerations, as is done in electromagnetism. In regard to this last point, damping 
proportional to velocity would not be expected to be physically relevant. While the 
dynamics ([!]) may be a reasonable description, for small q, of a free oscillator in a 
fluid, nano-oscillators with mechanical damping determined by the particular material 
geometry would not experience a simple damping proportional to velocity. Real 
macroscopic oscillators may however be describable by an effective susceptibility and 
by results such as the thermal energy (112). 
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Appendix A. Green function 



Here we present a Green function G(t,t ) that satisfies (32). We obtain the solution by 



solving (34) for G(u,to) with e > 0. The Green function corresponds to an infinitesimal 



value of e, but (34) is first solved for a non-infinitesmal number e > for which the 
standard solution method [31] can be applied without difficulty; after the solution is 
obtained the limit of infinitesimal e is taken and then G(u,to) is (inverse) Fourier 



transformed to find G(t, to). The solution of (34) with e > differs depending on whether 



to is positive or negative, as this affects the analytic properties of the inhomogeneous 



term as a function of co. The general solution of (34) will also contain the solution (25) 



of the homogeneous integral equation (24), which we do not require here and therefore 

and then take e — > + , 
0: 



drop. We find the solution G(u>, t ) by the standard technique 
whereupon the infinitesimal number + serves only to regulate a pole at co 



G(u,t ) 



(co + i 7 ) exp(io;tc 



(u - i0+)(u 2 -co 2 + iju) (u - i0+) [(co 2 -co 2 ) 2 + y 2 co 2 ] 



exp (-f ) [(co 2 - 7 >i cos (^) - 7 (7 2 



3coq) sin 



(T)] 



coi[(co 2 — col) 2 + 7 2 ^ 2 ] 



t >0, (A.1) 



G(u,t ) 



(co — iy) exp(io;to) 



[CO 



m)(co 2 



iyco) 



exp(f ) [(^- 7 > lC os(^)+ 7 (f 



3wq) sin 



m] 



loi[(lo 2 — col) 2 + 7 2 ^ 2 ] 



t < 0. (A.2) 



The solution (A.1)-(A.2) is continuous at to — 0, and co\ is again given by (27). Fourier 



transformation of (|A.1|)-(|A.2|) to the time domain gives the Green function 

Ji+t 



G(t,t ) = [28(-t) - 6(t 



t)) -V 

COn 



rxi)(-2M 



coicol 



7WiSgn(t) cos 



~2~ 



+ (7 2 - M 



sin 



LOit 

~2~ 



exp [- 



-0(t - to) 
exp [~l(\t\ + 1 



uicol 



2yuo 2 uo 2 



C0\ COS 



7CU1 COS 



LOit 

~2~ 



COi 



(t-to) +( 7 2 -2u; ( 



sin 



CO-i 



(t - 1 ) 



7 sin 



Ui\t\ 



G(t,t ) 



(y 2 — U)q)Ui cos 

= e(t 



Uit 



+ (7 2 -3w^)7 sin 



Uit 



tn > 0, 



(A.3) 



t)^e 0+ * 



0Oa 



-6{h - 1) 



exp [|(t — 1 ( 
coxcol 



yco\ COS 



(t-t ) -(Y-2co 2 )sm -^(t-t ) 



co x 
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(7 2 — Co>o Vi cos 



(7 2 - 3wq)7 sin 



to < 0, 



(A.4) 



where 9(x) is the step function. The values of G(t,to) at zeros of the step functions in 
(A.3)-( A~4| ) are to be taken as limits; G(t,t ) is in fact continuous across t — 0, to — 
and t = to. It is straightforward to verify that (A.3)-(A.4) satisfies (32); there is a 
discontinuity in dG(t, t )/dt at t = t which gives a delta function in d 2 G(t, t )/dt 2 . 



The infinitesimal number in (A.3)-(A.4) regularizes its Fourier transform in t, giving 



the G(u,t ) of (A.1)-(A.2). 



Appendix B. Coupling functions for a diagonalizable Hamiltonian 

Here we complete the final steps in diagonalizing the Hamiltonian. In section [5] 



we constructed a transformation, given by (61)-(63), (73), (75) and (79), that will 



diagonalize the Hamiltonian provided the commutation relations (52) and (53) are 



satisfied by (61)-(63). Inserting (61)-(63) into (52) and (53), and using (58), we obtain 



the final non-trivial conditions on the /-coefficients for the diagonalizing transformation: 



dw/gH/n (w) - / dw/!(w)/ n ,(w 



ih, 



du"f x {uj,uj")r nx {u'^")- 

du"f x (u>,u>")fx(u',u") - 
du"f nx (u,u")f^(u;\u") 





oo 



dco"f x (co,oj")f Ux ^'^") = ifoV 
du" f* x (u,u")f x {u' = 0, 
dJ'K{u,J')f Ux {u',J')=Q. 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



The last two of these conditions are easily seen to hold because of ( 73 ) and the first of 



(69). We proceed to show that the first two conditions, i.e. (B.l) and (B.2), hold for 



coupling functions satisfying (80) and (81), a class that includes the coupling function 



(21) corresponding to damping proportional to velocity. The following analysis has an 

although the results are different 



exact counterpart in the Huttner-Barnett model 
due to a different type of coupling between the g-oscillator and the reservoir 



Consider first the condition (B.l). If we insert the first of (68) together with (75) 
into (B.l), we encounter a term quadratic in h q (u). As noted after (76), a non-zero 



h q {u) will contain a delta function, so to avoid a square of a delta function in (B.l) we 
are forced to choose the solution 



h q (co) = 



of (76). The condition (B.l) is then, with use of (79) 
1 = 



dua (oj) 



u 2 - uj 2 Q + P 



, a 2 (Q ivra 2 M 



(B.5) 



(B.6) 



This complicated restriction on the coupling function a(u) is a requirement for the 
Hamiltonian to be diagonalizable. For a very broad class of coupling functions however, 
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(B.6) reduces to a much simpler requirement. If a 2 (u) is an even function of u then 



(B.6) can be written 
1 

1 = — 
in 



-lo 2 -P 



^ C2 



ma (tu 



£ 2 - u 2 2u 

The assumption that a 2 (u) is an even function also implies that the combination 



(B.7) 



+ 



17TQ! a; 



2w 



(B.8) 



which appears in the denominator in (B.7), is analytic on the upper-half complex em- 
plane. This follows from the familiar analysis of susceptibilities that are analytic on 
the upper-half frequency plane and that therefore obey Kramers-Kronig relations [33J; 
the real and imaginary parts of (B.8) exhibit one of the Kramers-Kronig relations for a 
complex susceptibility, provided a 2 (co) is an even function of u. If we further assume that 
a 2 (u)/uj > except possibly at u — 0, then it is shown in [33] that the "susceptibility" 
(B.8) does not take real values at any finite point in the upper-half complex cu-plane 
except on the imaginary axis, where it varies monotonically from its value at iO to its 
value at ioo. Moreover, the (real) value of (B.8) on the positive imaginary axis u = in, 
k > 0, is [33] 

« 2 (£) 







2 ' 



(B.9) 



We now consider the analytic properties of the integrand in (B.7) in the upper-half 
frequency plane. As (B.8) is real in the upper-half plane only on the positive imaginary 
axis, the denominator in (B.7) can vanish only at positive imaginary frequencies u = in, 
where it has the value (using (B.9)) 

« 2 (0 



U) Q + K 
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It is clear that (B.10) has a zero if and only if 

V(0 



Wn < 



(B.10) 



(B.H) 



If the equality holds in (B.ll) then (B.10) has a zero at lo = ik = 0, but this does not 
give a pole in the integrand in (B.7) unless (B.10) goes as n n , n > 1, as K — > 0. There 
are therefore no poles in the integrand in (B.7) in the upper-half plane if and only if 

« 2 (0 ( 

(B.12) 



either oj n > 



or uj. 
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and (|B.10h ~ K n , n < 1, as k ->■ 0. 



If (B.12) holds then the integral in (B.7) is equal to minus the integral over an infinite 
semi-circle in the upper-half plane. Taking u = Re 1 ^, R — )-oo,O<0<7ras the 
integration parameter along the infinite semi-circle, we find that the integral in (B.7) is 
iir, so that (B.7) is satisfied. 
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There remains the condition (B.2). We must substitute the first of (69), (73), 
(75), (B.5) and (79). This gives an integral relation that is again easy to evaluate for 
coupling functions for which a 2 (u) is an even function of u and a 2 (u)/u> > except 
possibly at uj = 0. Evaluation of the integral by closing the countour in the upper-half 
plane shows that the integral relation holds without restrictions on the coupling function 
beyond those already derived from (B.l). This analysis completes the proof that the 
Hamiltonian is diagonalizable for the class of coupling functions (80) if and only if (81) 
holds. 
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